This readme is is for the data submitted in conjunction with the paper
Title: "An ϵ-pseudoclassical model for quantum resonances in a cold dilute atomic gas periodically driven by finite-duration standing-wave laser pulses"
Authors: Benjamin T. Beswick, Hippolyte P. A. G Astier, Simon A. Gardiner, Ifan G. Hughes, Mikkel F. Andersen, and Boris Daszuta

These data were produced using the Python programing language and plotted using the Matplotlib.pyplot module. Figure 1 was produced using Inkscape.


========================================
Explanations of structure of data files
========================================


-------------------------
Files included in Figure 2/:
-------------------------

===============
(1) trajectory0.00
(2) trajectory0.05
(3) trajectory0.20
(4) trajectory0.25
===============

Each text file labled "trajectory*" (where "*" represents the value of the 𝛽 parameter for which the Poincare section was produced) contains a list of comma separated values of theta on the first line and a corresponding list of comma separated J values on the second. This data encompases all points plotted on each scatter plot including the initial phase space points. The data are ordered in such a way that the first set of (theta,J) values correspond to the first initial point plotted, and the subsequent values correspond to the trajectory resultant from this initial point, up to the total number of kicks plotted for each initial point. The values following this are the coordinates of the next initial point plotted, followed by those of its resultant trajectory, and so on.

-------------------------
Files included in Figure 3/:
-------------------------

===============
(1) trajectory0.25
(2) trajectory2.51
(3) trajectory5.01
(4) trajectory7.51
===============

Each text file labled "trajectory*" (where "*" represents the value of the ṽ parameter for which the Poincare section was produced) contains a list of comma separated values of theta on the first line and a corresponding list of comma separated J values on the second. This data encompases all points plotted on each scatter plot including the initial phase space points. The data are ordered in such a way that the first set of (theta,J) values correspond to the first initial point plotted, and the subsequent values correspond to the trajectory resultant from this initial point, up to the total number of kicks plotted for each initial point. The values following this are the coordinates of the next initial point plotted, followed by those of its resultant trajectory, and so on.

-------------------------
Files included in Figure 4/:
-------------------------

===============
(1) 0.1/momentumMomentPseudoclassical
(2) 0.1/momentumMomentQuantum
(3) 0.01/momentumMomentPseudoclassical
(4) 0.01/momentumMomentQuantum
(5) 0.63095734448/momentumMomentPseudoclassical
(6) 0.63095734448/momentumMomentQuantum
(7) 0.063095734448/momentumMomentPseudoclassical
(8) 0.063095734448/momentumMomentQuantum
(9) 0.158489319246/momentumMomentPseudoclassical
(10) 0.158489319246/momentumMomentQuantum
(11) 0.0158489319246/momentumMomentPseudoclassical
(12) 0.0158489319246/momentumMomentQuantum
(13) 0.251188643151/momentumMomentPseudoclassical
(14) 0.251188643151/momentumMomentQuantum
(15) 0.0251188643151/momentumMomentPseudoclassical
(16) 0.0251188643151/momentumMomentQuantum
(17) 0.398107170553/momentumMomentPseudoclassical
(18) 0.398107170553/momentumMomentQuantum
(19) 0.0398107170553/momentumMomentPseudoclassical
(20) 0.0398107170553/momentumMomentQuantum
===============

Each directory "*/", where "*" represents the value of the ϵ parameter for which the data were produced, contains two text files. Each file "momentumMomentPseudoclassical.txt" contains the line-break separated values of the pseudoclassical second order momentum moment, recorded after each kick for the value of ϵ given by the directory name. The files "momentumMomentQuantum.txt" contain the corresponding second order momentum moment calculated using the fully quantum treatment.


-------------------------
Files included in Figure 5/:
-------------------------

===============
(1) betaWaterfall0.001/0.0/momentumMomentPseudoclassical
(2) betaWaterfall0.001/0.0/momentumMomentQuantum
(3) betaWaterfall0.001/0.1/momentumMomentPseudoclassical
(4) betaWaterfall0.001/0.1/momentumMomentQuantum
(5) betaWaterfall0.001/0.2/momentumMomentPseudoclassical
(6) betaWaterfall0.001/0.2/momentumMomentQuantum
(7) betaWaterfall0.001/0.05/momentumMomentPseudoclassical
(8) betaWaterfall0.001/0.05/momentumMomentQuantum
(9) betaWaterfall0.001/0.15/momentumMomentPseudoclassical
(10) betaWaterfall0.2/0.15/momentumMomentQuantum
(11) betaWaterfall0.2/0.25/momentumMomentPseudoclassical
(12) betaWaterfall0.2/0.25/momentumMomentQuantum
(13) betaWaterfall0.2/0.0/momentumMomentPseudoclassical
(14) betaWaterfall0.2/0.0/momentumMomentQuantum
(15) betaWaterfall0.2/0.1/momentumMomentPseudoclassical
(16) betaWaterfall0.2/0.1/momentumMomentQuantum
(17) betaWaterfall0.2/0.2/momentumMomentPseudoclassical
(18) betaWaterfall0.2/0.2/momentumMomentQuantum
(19) betaWaterfall0.2/0.05/momentumMomentPseudoclassical
(20) betaWaterfall0.2/0.05/momentumMomentQuantum
(21) betaWaterfall0.2/0.15/momentumMomentPseudoclassical
(22) betaWaterfall0.2/0.15/momentumMomentQuantum
(23) betaWaterfall0.2/0.25/momentumMomentPseudoclassical
(24) betaWaterfall0.2/0.25/momentumMomentQuantum
===============

Each directory "betaWaterfall*/&", where "*" and "&" represent the values of the ϵ and 𝛽 parameter for which the data were produced, contains two text files. Each file "momentumMomentPseudoclassical.txt" contains the line-break separated values of the pseudoclassical second order momentum moment, recorded after each kick for the values of ϵ and 𝛽 given by the directory name. The files "momentumMomentQuantum.txt" contain the corresponding second order momentum moment calculated using the fully quantum treatment.

-------------------------
Files included in Figure 6/:
-------------------------

===============
(1) 0.2/populationMatrixQuantumZeroTemp
(2) 0.2/momentumMomentQuantumZeroTemp
(3) 0.2/populationMatrixClassicalZeroTemp
(4) 0.2/momentumMomentClassicalZeroTemp
(5) 0.2/populationMatrixQuantumFiniteTemp
(6) 0.2/momentumMomentQuantumFiniteTemp
(7) 0.2/populationMatrixClassicalFiniteTemp
(8) 0.2/momentumMomentClassicalFiniteTemp
(9) 0.02/populationMatrixQuantumZeroTemp
(10) 0.02/momentumMomentQuantumZeroTemp
(11) 0.02/populationMatrixClassicalZeroTemp
(12) 0.02/momentumMomentClassicalZeroTemp
(13) 0.02/populationMatrixQuantumFiniteTemp
(14) 0.02/momentumMomentQuantumFiniteTemp
(15) 0.02/populationMatrixClassicalFiniteTemp
(16) 0.02/momentumMomentClassicalFiniteTemp
(17) 0.11/populationMatrixQuantumZeroTemp
(18) 0.11/momentumMomentQuantumZeroTemp
(19) 0.11/populationMatrixClassicalZeroTemp
(20) 0.11/momentumMomentClassicalZeroTemp
(21) 0.11/populationMatrixQuantumFiniteTemp
(22) 0.11/momentumMomentQuantumFiniteTemp
(23) 0.11/populationMatrixClassicalFiniteTemp
(24) 0.11/momentumMomentClassicalFiniteTemp
===============

Each directory "*/", where "*" represents the value of the ϵ parameter for which the data were produced, contains eight text files. Each file "momentumMomentClassical^.txt" contains the line-break separated values of the pseudoclassical second order momentum moment, recorded after each kick for the value of ϵ given by the directory name. The files "momentumMomentQuantum^.txt" contain the corresponding second order momentum moment calculated using the fully quantum treatment. The files populationMatrixClassical^.txt and populationMatrixQuantum^.txt contain the full-sized arrays of momentum state populations calculated by our codes for the pseudoclassical and quantum models respectively, before the population cutoff was applied. The time axis rund from left to right in intervals of one kick, and each line break corresponds to a different momentum state, the central line is positioned at P=0 [hbar k]. The "^" symbol is replaced with "ZeroTemp" for the zero-temperature simulations and "FiniteTemp" for simulations run for 𝑤=2.5.

-------------------------
Files included in Figure 7/:
-------------------------

===============
(1) 0.2/populationMatrixQuantumZeroTemp
(2) 0.2/momentumMomentQuantumZeroTemp
(3) 0.2/populationMatrixClassicalZeroTemp
(4) 0.2/momentumMomentClassicalZeroTemp
(5) 0.2/populationMatrixQuantumFiniteTemp
(6) 0.2/momentumMomentQuantumFiniteTemp
(7) 0.2/populationMatrixClassicalFiniteTemp
(8) 0.2/momentumMomentClassicalFiniteTemp
(9) 0.02/populationMatrixQuantumZeroTemp
(10) 0.02/momentumMomentQuantumZeroTemp
(11) 0.02/populationMatrixClassicalZeroTemp
(12) 0.02/momentumMomentClassicalZeroTemp
(13) 0.02/populationMatrixQuantumFiniteTemp
(14) 0.02/momentumMomentQuantumFiniteTemp
(15) 0.02/populationMatrixClassicalFiniteTemp
(16) 0.02/momentumMomentClassicalFiniteTemp
(17) 0.11/populationMatrixQuantumZeroTemp
(18) 0.11/momentumMomentQuantumZeroTemp
(19) 0.11/populationMatrixClassicalZeroTemp
(20) 0.11/momentumMomentClassicalZeroTemp
(21) 0.11/populationMatrixQuantumFiniteTemp
(22) 0.11/momentumMomentQuantumFiniteTemp
(23) 0.11/populationMatrixClassicalFiniteTemp
(24) 0.11/momentumMomentClassicalFiniteTemp
===============

Each directory "*/", where "*" represents the value of the ϵ parameter for which the data were produced, contains eight text files. Each file "momentumMomentClassical^.txt" contains the line-break separated values of the pseudoclassical second order momentum moment, recorded after each kick for the value of ϵ given by the directory name. The files "momentumMomentQuantum^.txt" contain the corresponding second order momentum moment calculated using the fully quantum treatment. The files populationMatrixClassical^.txt and populationMatrixQuantum^.txt contain the full-sized arrays of momentum state populations calculated by our codes for the pseudoclassical and quantum models respectively, before the population cutoff was applied. The time axis rund from left to right in intervals of one kick, and each line break corresponds to a different momentum state, the central line is positioned at P=0 [hbar k]. The "^" symbol is replaced with "ZeroTemp" for the zero-temperature simulations and "FiniteTemp" for simulations run for 𝑤=2.5.

-------------------------
Files included in Figure 8/:
-------------------------

===============
(1) 0.1/momentumMomentPseudoclassical
(2) 0.63095734448/momentumMomentPseudoclassical
(3) 0.251188643151/momentumMomentPseudoclassical
(4) 1.58489319246/momentumMomentPseudoclassical
(5) 3.98107170553/momentumMomentPseudoclassical
(6) 10.0/momentumMomentPseudoclassical

===============

Each directory "*/", where "*" represents the value of the 𝑤 parameter for which the data were produced, contains a text file "momentumMomentPseudoclassical.txt", which contains the line-break separated values of the pseudoclassical second order momentum moment, recorded after each kick.

===================================================================
For queries, contact Ben Beswick at b.t.beswick@durham.ac.uk
